Proceedings of the DPF-2011 Conference, Providence, RI, August 8-13, 2011 



1 



Soft QCD measurements in the forward direction with the LHCb 
experiment 

A. Camboni, on behalf of the LHCb collaboration 

Departament de Estructura i Constituents de la l\/lateria, Universitat de Barcelona, Barcelona, Spain 



LHCb presents studies of particle production in minimum bias events in pp collisions at y's = 7 TeV. These 
studies include measurements of strangeness production, particle ratios, baryon-antibaryon ratios and charged 
particle production. The forward coverage and low pT acceptance of the experiment makes these measure- 
ments very complementary to those performed by the central detectors at the LHC. Further benefits arise from 
the powerful particle identification capabilities provided by the LHCb RICH system. The measurements are 
. . . compared with theoretical predictions. 
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CN ; 1. Introduction 

^: 

D ' LHCb is one of the four large experiments operating at the Large Hadron CoUider (LHC) at CERN. It is 
devoted to precision measurements of CP violation and rare B meson decays. Since the bb production in pp 
\^ ' collisions at 7 TeV is strongly favored in the forward/backward region, LHCb is build as a single arm magnetic 
C^l dipole spectrometer with a polar angular coverage with respect to the beam line of approximately 15 to 300 

mrad in the horizontal bending plane, and 15 to 250 mrad in the vertical non-bending plane. A full description 

^ of LHCb may be found in [1]. All subdetectors were fully operational and in a stable condition for the data 
that are analysed. For the measurements presented here the tracking detectors and particle identification (PID) 
JL , strategy are of particular importance. 

^ ■ The tracking system is composed by a precision Vertex Locator (VELO) surrounding the interaction region, 
' a dipole magnet and three downstream tracking stations. The resolution for primary (secondary) vertices is 
^ 50(100) /im and the momentum resolution for tracks having hits in both the VELO and the tracking stations 
^ ] is around 5%. Two Ring Imaging Cherenkov (RICH) detectors provide excellent charged particle identification 
^ I capabilities over a wide momentum range of 2 — 100 GcV/c. 

■ The results presented here are based on 6.8 /ib^^ from the LHC 2009 pilot run at y's = 0.9 TeV, 0.3 nb~^ 
, at ^/s = 0.9 TeV accumulated in early 2010 and 14 nb~^ accumulated by summer 2010 at y/s = 7 TeV. These 
^^O ■ data-sets were recorded using minimum bias triggers, which require a minimum energy deposit in the LHCb 
I calorimeters or at least one reconstructed track in the event. Due to the finite beam size and crossing angle, 
^\ the two halves of the VELO were partially retracted during the y/s = 0.9 TeV runs. 

, Soft QCD measurements are important input for modelling of the underlying event and tuning of event 
] ■ generators. Current models have been tuned to describe SPS and Tevatron data, with transverse momentum 
" . I I coverage of pt > 0.5 GeV and central rapidity. LHCb has the unique capability to explore uncovered phase 
^ . space regions. Good understanding of soft QCD processes is also required for extracting many important 
measurements at the LHC. 

The measurement presented in this paper are the production cross section of K*^ at -^/s = 0.9 TeV, the prompt 
d ' A/A, A/K^ and p/p production ratios, the inclusive (j> production cross section at i/i = 0.9 TeV and -^s — 7 
TeV and the measurement of prompt charged particle multiplicity. 
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2. Measurement of strangeness production 

Strangeness production studies provide sensitive tests of soft hadronic interactions, as the mass of the strange 
quark is of the order of the QCD scale Aqcd- Strange-hadron production is suppressed, as a consequence, but 
still occurs in the non-perturbative regime. Different fragmentation models, tuned on Tevatron data, agree on 
the total amount of strangeness produced but disagree on its distribution over phase space and on the ratios of 
the strange hadrons produced. 

LHCb is particularly interesting for QCD studies since one can measure the production of strange quarks in 
a range where models are expected to diverge more than for central rapidities. Moreover, V'^ decays {K^, A, 
A) have a very clean experimental signature which allows to identify them unambiguously using only kinematic 
information. As for the (j) meson, it offers a unique way to study primary strangeness production, and via its 
decay into charged kaon pairs is vital for the calibration of the PID system of LHCb. 
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2.1. K° production cross section 



The hadronic production of mesons has been studied by several experiments at a range of different centre- 
of-mass energies, both in pp and pp collisions (see for example 0, Q)- The most recent measurements of 
production at the Tevatron have shown deviations with respect to the expectations of hadronization models 0] . 

This is an ideal early measurement since K'^ are relatively abundant in minimum bias data and K'^ — > tt+tt^ 
candidates from the 2009 pilot run can be selected using only the tracking system. A high-purity sample could 
then be obtained without using particle identification information. 

In Figure [U the LHCb results is presented in bins of y and pT- The different Monte Carlo (MC) predic- 
tions agree reasonably well with the data, although they tend to underestimate (overestimate) the measured 
production in the highest (lowest) px bins. 
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Figure 1: Double-differential prompt K'^ production cross-section in pp collisions at -^s — 0.9 TeV as a function of 
transverse momentum pt and rapidity y. The points represent LHCb data. Measurements are compared with the 
Pythia tune PerugiaO |4|] and the LHCb MC with and without diffractive events. The lower plots show the MC/data 
ratios, with the shaded band representing the uncertainty for one of these ratios. 



Due to the long lifetime and partially open VELO position, only a small fraction of the daughter 
tracks traversing the spectrometer leave a signal in the VELO. Therefore, a twofold path was followed for the 
reconstruction and selection, on one hand using tracks reconstructed only with hits in the tracking stations, 
on the other hand forming candidates with tracks leaving hits also in the VELO. The results are then 
combined. All the efficiencies are estimated using MC simulations. 

To determine the absolute luminosity, for this study, a novel method was designed 0, Q , employing the high 
resolution of the LHCb vertex detector to directly measure width and position of the beams using beam-gas 
and beam-beam interactions. Information about the bunch currents was provided by the LHC machine. In 
this analysis the main systematic contributions were given by the bunch intensities, 12%, and by the tracking 
efficiency, 10% in the most challenging bins. 

The LHCb measurement is the first one at 0.9 TeV. As one can see from Fig. [21 compared to earlier 
measurements the range in rapidity, y, and transverse momentum, px , were extended towards larger and 
smaller values, respectively. In general the agreement with previous measurements is reasonable, given the 
range of center-of-mass energies and that the results are averaged over different ranges in rapidity or transverse 
momentum. A detailed description of the analysis procedure and results are reported in fsl . 



Proceedings of the DPF-2011 Conference, Providence, RI, August 8-13, 2011 



3 



^ 10 

s 

E 



Q. 



0.1 



0.0 



4+ 





LHCb pp 


900 GeV 


2.5 < y < 4.0 





UAl p p 


630 GeV 


-2.5 < 11 < 2.5 


T 


CDF p p 


630 GeV 


-1 .0 < y < 1 .0 


A 


CDF p p 


1800 GeV 


-1.0< y< 1.0 


1 


UA5 p p 

, , 1 , 


540 GeV 
, , 1 , 


-3.5 < y<3.5 
, , 1 , 


0.2 


0.4 


0.6 


0.8 



1.0 



1.2 



. I 
1.4 



III. 
1.6 

[GeV/c] 



Figure 2: Absolute measurements of the Ks production cross-section as a function of pr, performed by tiie UAl, UA5, 
CDF and LHCb experiments, at different liigli-energy liadron colliders and in different y or r] ranges. 



2.2. Inclusive cross section 

A measurement of the inclusive differential cj) cross-section in pp collisions at -y/s = 7 TeV is presented in this 
study mesons are reconstructed using the K^K~ decay mode and thus the selection relies strongly on 
LHCb's RICH detectors for particle identification purposes. 

Two data samples are used, one in which at least one kaon is required to pass a RICH based selection, the tag 
sample, and one in which both kaons have to pass this selection, the probe sample. To minimize the dependence 
on MC simulations the PID efficiency in , j/ bins is estimated with this tag and probe technique. The LHCb 
luminosity used in this cross-section estimation is based on a continuous analysis of hits in the Scintillator Pad 
Detector which has been normalized to the absolute luminosity scale. The absolute scale was determined 
with the beam-gas method mentioned before. The probability for multiple pp collisions per bunch-crossing was 
negligibly low in the data taking period considered. 

In Figure |3] are shown projections on the y and pt axis within the same kinematic region. Also plotted are 
the results from the Perugia 3 tuning of Pythia 6.4 and the standard LHCb MC [s!] based also on the same 
version of Pythia. In general, both Monte Carlo samples underestimate the measured meson production in the 
phase space region of the measurement. There is a tendency to better agreement for lower transverse momenta. 
The cross-section for inclusive production in the kinematic range 0.6 < pr < 5.0 GeV/c and 2.44 < y < 4.06 
is found to be: 

a{pp (1>X) = (1758± 19(stot) tfi{syst) ± im{scale)) fih . 



3. Baryon number transport and baryon suppression 

Other probes of strangeness production in high energy hadron collisions are cross-section ratios, where lu- 
minosity and many systematic uncertainties cancel. In baryon- antibaryon ratios like A/A or p/p, the baryon 
contains valence quarks in common with the proton, while all three antiquarks of the antibaryon have to be 
produced in the collision. The ratio of the production cross sections thus measures the baryon-number transport 
from the beam particles to the final state. Several models exist to describe this transport, but the mechanism 
is not well understood. Also the strange A/K^ ratio is a good test for different fragmentation models, which 
do not agree between them in this phase space region (9|] and where there were no experimental data available 
before. 

For the V'^ studies, high-purity prompt K^J , A and A samples were selected based on a Fisher discriminant 
constructed using the logarithms of the impact parameters of the particles and of their daughters ■ The 
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Figure 3: Inclusive differential <f) production cross-section as a function of pt (left) and y (rigfit), measured with data 
(black), and compared to the LHCb MC (solid, red) and Perugia (dashed, blue). The inner error bars represent the 
statistical uncertainty, the outer show the quadratic sum of statistical and systematic errors. Error bars shown in the 
ratio plots are statistical only. 



selection requirement of a primary vertex ensured that mainly the candidates coming from non-diffractive events 
are kept. For the two available energies the analyzed data correspond to 0.31 nb~^ at yfs = 0.9 TeV and 1.8 
nb"^ at = 7 TeV. 

To compare results at both collision energies, and to probe scaling violation, both production ratios are shown 
in Fig. |4]as a function of rapidity loss. Ay — j/bcam — y where ybcam is the rapidity of the protons in the anti- 
clockwise LHC beam, which travels along the positive z direction through the detector. Excellent agreement is 
observed between results at both -/s — 0.9 and 7 TeV as well as with results from STAR at ^/s = 0.2 TeV. The 
measured ratios are also consistent with results published by ALICE and CMS. 
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Figure 4: The ratios (a) A/A and (b) K/K° from LHCb are compared at both t/s = 0.9 TeV (triangles) and 7 TeV 
(circles) with the published results from STAR [3| (squares) as a function of rapidity loss. Ay — j/beam — y- Vertical 
lines show the combined statistical and systematic uncertainties and the short horizontal bars (where visible) show the 
statistical component. 



Both measured ratios are compared to the predictions of the Pythia 6 generator tunes: LHCb MC, Perugia 
and Perugia NOCR ^J], as functions oipr and y at y/s = 0.9 and 7 TeV (Fig. [5]shows the latter). According 
to Monte Carlo studies, the requirement for a reconstructed primary vertex results in only a small contribution 
from diffractive events to the selected V° sample, therefore non-diffractive simulated events are used for these 
comparisons. The predictions of LHCb MC and Perugia are similar throughout. The ratio A/A is close to 
Perugia at low y but becomes smaller with higher rapidity, approaching Perugia NOCR. In collisions at -/s 
— 7 TeV, this ratio is consistent with Perugia across the measured pr range but is closer to Perugia NOCR 
at = 0.9 TeV. The production ratio A//f° is larger in data than predicted by Perugia at both collision 
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energies and in all measurement bins, with the most significant differences observed at high pT ■ 

A slight increase in the ratio ^/K^ is observed when going from y/s = 0.9 TeV to 7 TeV. This is plausible, since 
particle masses and kinematic factors in general should become less important at higher energies. Since both 
the A and the contain a single strange valence quark, the discrepancy cannot be explained by a mismatch in 
strangeness suppression between data and Monte Carlo. Instead it probes directly the understanding of baryon 
production in the fragmentation process. 
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Figure 5: The ratios A/A and A/K^ at -^5 = 7TeV compared with the predictions of the LHCbMC, PerugiaO and 
Perugia NOCR as a function of (a) & (c) rapidity and (b) & (d) transverse momentum. Vertical lines show the combined 



statistical and systematic uncertainties and the short horizontal bars (where visible) show the statistical component. 



For the measurement oip/p ratio a data sample similar to the previous study was used RICH information 
allows to isolate almost pure samples of pions, kaons and protons using delta log-likelihood cuts on particle type 
hypotheses. The hadron particle identification was calibrated using secondaries from the decays — >■ 7r"'"7r~ 
and A — > ptt~ , with pions and protons identified using kinematic cuts only. Charged kaons were selected from 
the decay (j) — ^ K~ using a tag and probe approach where only one of the kaons was required to be identified 
by RICH. Similar to the baryon transport pattern for A/A, the p/p ratio shows a big deviation in ratio from 
unity at low energy (except for the highest pr bin). Fig. |6l These results are compared against the predictions 
of Pythia as implemented in the LHCb simulation and those of the so-called Perugia tune of Pythia [3| . In this 
study LHCb Pythia includes the simulation of diffractive events, whereas Perugia does not. For -^i = 0.9 TeV 
at high pt there is good agreement between the data and the Perugia predictions, but less so for LHCb Pythia; 
in the low and middle px bins at this energy the data values are in general below those of both generators. At 
= 7 TcV the measurements are compatible with the expectation of both LHCb Pythia and Perugia 0. The 
results at both energies can also be considered in terms of the rapidity loss Ay = yhca.ui — y, the measurements 
are shown in this form in Fig. [T] These span an interval in Ay of almost 4 units, reaching up to a highest 
value of close to Ay = 7. In a given pT bin the measurements appear consistent with a monotonic distribution. 
There is an indication that the results have some px dependence. Figure [7] also shows measurements of the 
same quantity performed by other experiments. Reasonable consistency is observed. 
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Figure 6: Distribution of the ratio p/p against y: up at ^/s — 0.9 TeV, down at ^/s — 7 TeV. 
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Figure 7: LHCb measurements of the p/p ratio plotted versus Ay. 



4. Measurement of Charged Particle Multiplicities 

The charged particle production in proton-proton collisions at a centre-of-mass energy of ^/s — 7 TeV in 
different bins of pseudo-rapidity rj have been studied with the LHCb detector and compared with model predic- 
tions. The charged particles are reconstructed in the VELO detector close to the interaction region, essentially 
outside of the magnetic field, which provides high reconstruction efHciency in the rj ranges —2.5 < rj < —2.0 
and 2.0 < rj < 4.5. In the absence of almost any magnetic field in the VELO region, no further acceptance 
corrections as function of momentum are needed. Since the VELO detector is close to the interaction region, the 
amount of material before the particle detection is minimal, minimizing the corrections for particle interactions 
with detector material. The data were taken with a minimum bias trigger, only requiring a minimum of one 
reconstructed track in the VELO detector. 
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Monte Carlo event simulation is used to correct for acceptance and resolution effects. The detector simulation 
is based on the Geant 4 [l^l program and incorporates our best knowledge of the LHCb spectrometer. The 
Monte Carlo event samples were passed through reconstruction and selection procedures identical to those for 
the data. Proton-proton collisions, both elastic and inelastic processes, were generated using the Pythia 6.4 
event generators which was tuned to lower energy hadron collider data [8| . 

The reconstructed multiplicity distributions are corrected on an event by event basis to account for the 
tracking and selection efficiencies and for the background contributions. These corrected distributions are then 
used to measure the charged particles multiplicities in each of the rj ranges though an unfolding procedure. The 
distributions are corrected for pile- up event. Hard interactions can be enriched by requiring at least one high 
Pt > IGeV/c track in the forward acceptance of the event. The acceptance is no longer independent of the 
momentum and requires an additional correction to take into account the efficiency to reconstruct a high px 
track. 

The LHCb Monte Carlo is used to estimate the overall tracking and selection efficiency as a function of 
pseudorapidity and azimuthal angle (j). The main correction come from non-prompt particle contamination 
(5-10%), mainly tracks from converted photons. 

Fig. [S] shows the charged particle multiplicity distribution for different ranges in pseudorapidity rj. In general 
all generators underestimate the multiplicity distributions, with the LHCb tune giving the best description of 
the data. The exclusion of the Pythia diffractive processes in the Perugia tunes also improves the description 
of the data, particularly in the full forward region. 

The charged particle density as a function of pseudorapidity is shown in Fig. [5] The data has a marked 
asymmetry between the forward and backward region; this is a consequence of the requirement of at least one 
track in the full forward range. All models fail to describe the mean charge particle multiplicity per unit of 
pseudorapidity. 
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5. Conclusions 



The first studies with data from the 2009 and 2010 runs showed that the LHCb experiment is ready for its core 
physics program. Early data were used to perform and to check the calibration of the different subdetectors, 

as well as to perform new measurements in a unique range of rapidity and transverse momentum. A first 
cross-section measurement at 0.9 TeV was produced at larger rapidities and smaller transverse momenta than 
previous measurements. Preliminary results for rj inclusive cross section were shown. Preliminary measurements 
for ratios of V'^ and protons suggest lower baryon suppression and higher baryon transport in data than in the 
Monte Carlo models investigated. The LHCb data are consistent with data from lower energy experiments. In 
genciral, the models underestimate the charged particle production, LHCb has already provided great input for 
their improvement. 
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